Three-body correlations in Borromean halo nuclei 
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Three-body correlations in the dissociation of two-neutron 
halo nuclei are explored using a technique based on intensity 
interferometry and Dalitz plots. This provides for the com- 
bined treatment of both the n-n and core-n interactions in 
the exit channel. As an example, the breakup of 14 Be into 
12 Be+n+n by Pb and C targets has been analysed and the 
halo n-n separation extracted. A finite delay between the 
emission of the neutrons in the reaction on the C target was 
observed and is attributed to 13 Be resonances populated in 
sequential breakup. 



The quest for the drip- lines, which define the limits of 
binding for nuclear systems, has only been attained for 
light nuclei. As such, these nuclei are unique in displaying 
the manner in which nucleons bind in an A nucleon sys- 
tem from the most neutron deficient to the most neutron 
rich. Clustering phenomena, observed for example re- 
cently in excited states close to threshold [Q , also appear 
as haloes in ground states near the neutron drip-line Q . 
The most intriguing manifestation of clustering are the 
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Borromean two-neutron halo nuclei ( 6 He, n Li and 14 Be), 
in which the two-body subsystems are unbound || . Be- 
yond the global properties, such as the abnormal sizes 
[4,|| or the low momentum content of the constituents 
[5|, the experimental challenge lies in determining the 
structure of these three-body systems. 

The dissociation in the field of a target nucleus, fol- 
lowed by the measurement of the momenta of the frag- 
ments (core+n+n), has been used in attempts to probe 
correlations in two-neutron halo nuclei. Vestiges, how- 
ever, of the two-body forces that stabilise the projectile 
in the ground state may affect the three-particle decay 
in the form of final-state interactions (FSI). Experimen- 
tally, beyond the reconstruction of the core+2n invariant 
mass, the analyses so far reported have been restricted to 
the binary channels. The relative energy in the core+n 
channel has been used, for example, to probe the for- 
mation of core-n resonances The n-n observables, 
which should provide access to correlations within the 
halo, have often been compared only to simplified inter- 
pretations, such as a di- neutron configuration or three- 
body phase space neglecting the n-n FSI. 

In the present paper the three-body correlations in the 
dissociation of two-neutron halo nuclei are explored. In 
particular, a new method for analysing triple coincidence 
events (core+n+n) from kincmatically complete experi- 
ments is described. The method incorporates the tech- 
niques of intensity interferometry [ [TT| and Dalitz plots 
|12| and permits the halo n-n separation and time delay 
between the emission of the two neutrons to be derived. 
As will be seen, the latter is related to the presence of 
core-n FSI in the exit channel. In principle, the present 
approach also allows the energies and lifetimes of these 
resonances (or virtual states if I = 0) to be derived. 

When neutrons are emitted in close proximity in space- 
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time, the wave function of relative motion is modifie d by 
the known FSI and quantum statistical symmetries . 
Two-neutron intensity interferometry, and in particular 
the correlation function | ft3f C nn , relates this modification 
to the space-time separation of the particles at emission: 



Cnn(pi,P2) 



d 2 n/ dp\dp2 
(dn/dpi) (dn/dpa) 
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The numerator is the measured two-particle distribution 
and the denominator the product of the independent 
single-particle distributions |l3| , |l4| , which are commonly 
projected onto one dimension, the relative momentum 
q = \pi — p2 1 . In an earlier paper (l4| , the correlation 
function was extracted for the dissociation of 14 Be by Pb 
with the aim of probing the spatial configuration of the 
halo neutrons. The analysis followed the formalism of 
Rcf. |u| and assumed that the neutrons were emitted si- 
multaneously following dissociation in the Coulomb field 
of the target. A n-n separation of r™ s = 5.6 ± 1.0 fm 
was thus obtained. 

The same analysis has been applied to dissociation of 
14 Be by a C target, in order to investigate the influence 
of the reaction mechanism Jl6| ] . The n-n separation ob- 
tained was somewhat larger, r™ s = 7.6 ± 1.7 fm. This 
raises the question as to whether simultaneous emission 
can be assumed a priori. In principle, the analysis of the 
correlation function in two dimensions, transverse and 
parallel to the total momentum of the pair, would al- 
low for the unfolding of the source size and lifetime [^5) . 
Such an analysis requires a large data set and was thus 
not applicable to the present measurements. The sys- 
tem being studied here is far less complex, however, than 
those usually encountered in interferometry (for example 
compound nuclei evaporating particles or systems of col- 
liding heavy ions |ll]]). Moreover, the simple three-body 
nature of the system breaking up suggests immediately 
that any delay in the emission of one of the neutrons 
will arise from core-n FSI/resonances in the exit chan- 
nel (Fig. 0). As will be seen later, the degree to which 
such resonances are present may depend on the reaction 
mechanism. 
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FIG. 1. Schematic view of the sequential breakup of 
a two-neutron halo nucleus, whereby dissociation passes 
through a resonance in the core-n system. A delay r nn , equiv- 
alent to the lifetime of the resonance, is introduced in the 
emission of the second neutron. The effective source size is a 
convolution of the n-n separation r nn and the delay r nn . 



The data examined here were acquired from the dis- 
sociation of a 35 MeV/N 14 Bc beam into 12 Be+n+n by 
Pb and C targets. Details of the experiment and previ- 
ous analyses have been reported elsewhere 0,0. Each 
event was reconstructed from the momenta pjo of the 
neutrons measured using the DEMON array [|l8| as fol- 
lows: (i) we calculate the average velocity ((3) of the 
core+n+n frame at dissociation as that for which the 
mean for all events of the total neutron momentum along 
the beam axis (pi +j>2)z is and the average decay energy 
of the system (E^) (see below) is a minimum; (ii) in this 
frame, momentum conservation is applied event-by-event 
to reconstruct the core momentum p*3 = — (pi+p*2). From 
the four-momenta of the three particles wc calculate the 
total energy available in the center-of-mass of the system 
and extract the kinetic part, the decay energy: 



Ea = 
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Note that in the present analysis the decay has been re- 
constructed from only the neutron momenta pi, 2' The 
use of this reconstruction algorithm instead of the more 
classical analysis employing the core momentum p*3 mea- 
sured in the charged-particle telescope eliminates the lim- 
ited energy and position resolution of this detector and 
the uncertainty related to the depth within the target at 
which the reaction took place. Importantly, the form of 
the spectra obtained here for the 14 Be and core-n invari- 
ant masses agree well with those obtained using the core 
momentum p"7] , |l9| ]. 

An interacting phase-space model has been developed 
for the analysis of triple correlations in the data. In brief, 
the experimental decay energy distribution is used as in- 
put to generate events pi,2,3(-Ed) following three-body 
phase space pfjfl . The core-n resonances are introduced 
following the sequential breakup of the system (Fig. |l|) 
into one neutron and the core-n resonance with a relative 
energy E 2 z given by a Breit-Wigner distribution (E ,F); 
the resonance is then allowed to decay into the core plus 
neutron. In the n-n channel, the FSI is introduced via 
a probability P(\pi — p%\) to accept the event following 
the form of the measured n-n correlation function [0-|l6j . 
The final momenta pi, 2, 3 are filtered through a simulation 
including all experimental effects [QT 18 and the recon- 



struction algorithm described above is applied to pi ,2- 

Correlations in three-particle decays have been exten- 
sively studied in particle physics by means of Dalitz plots 
of the particle energies (Ei , Ej ) or the squared invariant 
masses of particle pairs (M?j, M? k ), with Mfj = (pi+pj) 2 . 
In these representations, FSI/resonances lead to a non- 
uniform population of the surface within the kinematic 
boundary defined by energy-momentum conservation and 
the decay energy [Q. The classic example of such an 
analysis is the three-body decay of an unstable particle 
pl[ . In the present case, the core+n+n system exhibits 
a distribution of decay energies. Consequently, the value 
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of Ed associated with each event will lead to a different 
boundary for the Dalitz plot, and the resulting plot con- 
taining all events cannot be easily interpreted. We have 
thus introduced a normalized invariant mass: 



m : 



Mf 3 



(rrii + rrij) 2 
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which ranges between and 1 (E^j between and Ed) 
for all events and exhibits a single kinematic boundary. 
Examples of how n-n and core-n FS1 present in the de- 
cay are manifested in core-n versus n-n Dalitz plots are 
displayed in Fig. ^ whereby events have been simulated 
with the model described above. The inputs were an Ed 
distribution following that measured (Fig. ||), the C nn 
obtained with the C target (Fig. ||), and core-n res- 
onances with T = 0.3 MeV at E = 0.8, 2.0, 3.5 MeV. 
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FIG. 2. Dalitz plots (core-n versus n-n) for simulations 
with the interacting phase-space model of 14 Be dissociation: 
without FSI (a), with n-n FSI (b) and with a core-n resonance 
at E = 0.8 (c), 2.0 (d) and 3.5 MeV (e). The combination of 
the n-n and the core-n FSI of (c) is shown in (f). 

In the absence of any FSI, the Dalitz plot exhibits, 
as noted above, a uniform population (Fig. |2p). The n-n 
FSI appears as a concentration of events with m^ n < 0.25 
)), which correspond to small relative momenta 
The core-n resonance at E — 0.8 MeV (Fig. ||c) 
appears as horizontal bands around m^ n ~ 0.25 and 0.75. 
The location of these bands depends on the energy of the 
resonance with respect to the mean decay energy of the 
system: a single band at m^ n s» 0.5 if Eq ~ (Ed) (Fig. Hd ) 
and two symmetric bands if E a ^ (Ed) (Fig. |p,e) fllq . 
This feature arises as when one of the neutrons forms 
a resonance with the core at a given value of m^, the 
relative core-n energy of the other is essentially fixed at 




m?, « 1 — |16|. Importantly, the Dalitz plot repre- 



13 ~ J- '"-23 

sentation provides not only for the identification of the 
different FSI between the particle pairs, but also a di- 



rect comparison of the relative importance of each in the 
decay (Fig. |f). 
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FIG. 3. Dalitz plots (core-n versus n-n), and the projec- 
tions onto both axes, for the data from the dissociation of 
14 Be by Pb (upper) and C (lower panels). The lines are the 
results of the phase-space model simulations with/ without 
(solid/dashed) n-n FSI. The insets correspond to the 14 Be 
decay energy. 

The Dalitz plot for the data from the dissociation by 
Pb (Fig. @) presents a strong n-n FSI and a uniform den- 
sity for m^ n > 0.5. Indeed, the n-n FSI alone describes 
very well the projections onto both axes, and therefore 
suggests that core-n resonances are not present to any 
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significant extent. This result confirms the hypothesis of 
simultaneous n-n emission employed in the original anal- 
ysis of the dissociation of 14 Be by Pb flJJ] . The value of 
r™ s so extracted, 5.6 ±1.0 fm, thus corresponds to the 
n-n separation in the halo of 14 Be. 

For dissociation by the C target (Fig. |^), despite the 
lower statistics, two differences are evident. Firstly, the 
n-n signal is weaker, indicating, as discussed earlier, that 
a significant delay has occurred between the emission of 
each neutron. Second, and more importantly, the agree- 
ment between the model including only the n-n FSI and 
the data for m^ n is rather poor. In order to verify whether 
this disagreement corresponds to the presence of core-n 
resonances, which would be responsible for the weakening 
of the n-n signal, we have investigated the core-n relative 
energy, E cn . It has been reconstructed for the simula- 
tions incorporating only the n-n FSI and compared in 
Fig. H to the data (the model calculations have been nor- 
malized to the data above 4 MeV). For dissociation by 
Pb, the inclusion of only the n-n FSI provides a very 
good description of the data, with the exception of small 
deviations below 1 MeV. This is in line with the Dalitz 
plot analysis discussed above. 

The deviations observed for the C target between the 
measured m^ n and the simulation including only the n-n 
FSI (Fig. H) clearly correspond to structures in the -B C n 
spectrum. Moreover, these structures are located at en- 
ergies that are in line with those of states previously re- 
ported in 13 Be: the supposed <fc/2 resonance at 2.0 MeV 
jp2f and a lower-lying state(s) |1^ , |22],^3|| , also suggested 
by various theoretical calculations ]24|-|27|]. The model- 
to-data ratio is about 1/2, indicating that the peaks cor- 
respond to resonances formed by one of the neutrons in 
almost all decays; the solid line in Fig. ^ accounts for the 
contribution of the neutron not interacting with the core. 
If we add to the phase-space model with n-n FSI core-n 
resonances (r = 0.3 MeV) for all events at E = 0.8, 2.0 
|| and 3.5 Me\|] with intensities of 45, 35 and 20%, 
respectively, the data are well reproduced (dashed line). 
In the case of dissociation by Pb, the lowest-lying reso- 
nance^) appears to be present in at most 10% of events. 

The different results obtained for the Pb and C targets 
may be attributed to the associated reaction mechanisms 
p8[ . In the case of the Pb target , the dominant process is 
electromagnetic dissociation |L7|], whereby the halo neu- 
trons behave as spectators and only the charged core is 
acted on by the Coulomb field of the target [|9| . Quali- 
tatively then, the n-n FSI may be expected to influence 
most strongly the decay. In the case of the C target, nu- 
clear breakup dominates and the reaction takes place at 



x The present data are not particularly sensitive to the lo- 
cation and form of the states, in particular below 1 MeV, 
and a resonance at 0.5 MeV would, for example, equally well 
describe the data. 



smaller impact parameters, in general through the inter- 
action of one of the halo neutrons with the target |fl7| . 
As such the population of core-n resonances is favoured 
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FIG. 4. Core-n relative energy distributions (upper) and 
n-n correlation functions (lower panels) for the dissociation of 
14 Be by Pb and C. The lines in the E cn spectra are the result 
of the phase-space model simulations with n-n FSI (solid) 
plus core-n resonances (dashed, see text). The histograms 
presented in the middle panels are the difference between the 
data and the n-n FSI simulations. The solid lines in the lower 
panels are the C nn for r™ s = 5.6 fm and r nn = 0; the shaded 
area corresponds to the inclusion of a finite r nn (see text). 

By combining the information extracted from the core- 
n channel with the n-n correlation functions, we can ex- 
tend the analysis and also extract the average lifetime of 
the core-n resonances. If we fix the n-n distance in 14 Be 
as that obtained for dissociation by Pb, r™ s = 5.6 fm, we 
can introduce the delay between the emission of the neu- 
trons T nn needed to describe the n-n correlation function 
for the C target. As discussed earlier, this delay should 
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correspond to the lifetime of the resonances (Fig. |l|) . The 
result, shown as the shaded area in Fig. f|, suggests an 
average lifetime of 150±i°q W c > or x 10 ~ 22 s - 

In summary, three-body correlations in the dissocia- 
tion of two-neutron halo nuclei have been explored. A 
new analysis technique employing intensity interferome- 
try and Dalitz plots has been presented and applied to 
the breakup of 14 Be by Pb and C targets. Through the 
combined treatment of both the n-n and core-n corre- 
lations, the halo n-n separation has been extracted and 
a finite delay found between the emission of the neu- 
trons for the reaction on C. This delay can be attributed 
to resonances in 13 Be populated in sequential breakup. 
The application of the techniques presented here to a 
well established system such as 6 He would be of particu- 
lar interest, as would the investigation of multi- neutron 
haloes. Finally, as the technique of intensity interferome- 
try is applicable to protons, proton-rich nuclei exhibiting 
similar few-body clustering may also be explored. 
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